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ABSTRACT 


Motivation to expand the understanding of a helicopter 
rotor descending into the vortex ring state (VRS) stems 
from the aircraft mishaps that have plagued the helicopter 
community. The V-22 has become the most recent victim of 
encounters with VRS. 

The onset of VRS is associated with the collapse of 
the helical vortex wake in the plane of the rotor. The 
resulting wake disturbances develop an irregular and 
aperiodic flow. Rotor blade interaction with the disturbed 
vortices causes large variations in the blade spanwise 
aerodynamic load distribution. Harmonic analysis of the 
loading indicates that higher harmonic content becomes 
prevalent in this state. 

The dynamic flow similarities achieved in a water 
tunnel are used to explore flow visualization and conduct 
vibration analysis of a rotor system operating in the VRS. 
A scaled rotor system was operated in the NFS Aeronautical 
Engineering Department's water tunnel. Sensors were used to 
gather thrust and vibration power spectrum data when 
operating in VRS. Experimental results correlate with full 
scale flight data and show a significant increase in the 
vibration levels of the even multiples of the blade passage 
frequency. The relative strength of these higher harmonics 
can be used as an indicator of impending VRS encounters. 
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INTRODUCTION 


I. 

A. MOTIVATION FOR RESEARCH 

The need for increased understanding of the 
aerodynamic effects of the vortex ring state (VRS) has 
become an issue of great importance, since the mishap 
involving the tilt rotor V-22 Osprey and the primary causal 
factor determined as asymmetrical vortex ring state. In 
fact, thirty-four documented helicopter incidents involving 
encounters with the vortex ring state (VRS) occurred 
between 1982 and 2003. [Varnes, 1999] 

To expand the understanding of this flow phenomenon, 
the development of a hydrodynamic facility to explore the 
flow field of a rotor system operating in the vortex ring 
state was undertaken. An experiment was conducted in the 
Naval Postgraduate School Aeronautical Engineering 
Department water tunnel, taking advantage of flow 
similarity factors, which permitted reduced operating 
velocity. The limited research conducted using a water 
tunnel also motivated studying this flow phenomena in a new 
medium. This thesis demonstrates that VRS boundary limits 
can be mapped and relevant vibration analysis can be 
performed on rotor systems operating in a water tunnel. 

B. HELICOPTER ROTOR FLOW STATES 

A rotor in vertical flight is typically described as 
being in one of three flow states. These are: (1) propeller 
working state; (2) vortex ring state; and (3) windmill 
brake state. 

1. Propeller Working State 

The propeller working state is terminology used to 
indicate vertical climb conditions. Figure 1. Airflow is 
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through the rotor in the opposite direction of thrust. 
During the propeller working state the trailing tip 
vortices are downstream of the path of each subsequent 
passing rotor blade. The limiting condition of the 
propeller working state is the hover. [Stewart, 1951] 



Figure 1. Helicopter Flow States. [From: Prouty, 

1986] 

2. Vortex Ring State 

The blade tip vortices, which normally cause minor 
losses in rotor efficiency, can be drastically enlarged 
when the aircraft begins to descend into its own downwash. 
Figure 1. Blades begin to recirculate the downwash quickly 
back above the rotor disk, only to be recirculated back 
down again. The airflow through the rotor plane does not 
have a regular slipstream, but flow is thrust out radially 
and recirculated. The rotor blades can eventually begin to 
stall as the angle of attack is increased due to the upward 
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airflow. Prior to blade stall, the recirculated air can 
almost be considered a closed system. Blade tip vortices 
from the rotor disk become trapped in this recirulation 
ring shape enclosing the outer rim of the rotor disk. If 
additional blade pitch angle is added, corresponding to an 
increase in collective, the rotor only recirculates the 
closed system blade wash at an increasing rate. 

Theoretically, the vortex ring state begins at rates 
of descent just below hover and transitions into the 
windmill brake state. From a practical standpoint, extreme 
flow variations do not start until the rate of descent is 
about half the hover induced velocity. Hover induced 
velocity is that downwash velocity which occurs when rotor 
thrust equals the gross weight of the aircraft. 

A well defined slipstream ceases to exist as an 
aircraft begins to descend. The region from descent to the 
windmill brake state can be referred to as the vortex ring 
state. A more commonly accepted value at which onset of the 
vortex ring state occurs is when the rate of descent equals 
0.7 times the induced velocity at hover. 
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Figure 2. Recirulation During Flight in the VRS. 

[From: Drees and Hendal Film, 1951] 

The vortex ring state has been of concern to the 
designers of helicopters since the 1950's. This phenomenon 
was first recognized by George de Bothezat in 1922. Early 
studies to verify the existence and its effects yielded 
results that are consistent with today's data. 

Rotor operations in the vortex ring state will cause 
some of the most violent aircraft vibrations of any flight 
regime. Flight conditions are characterized by high rates 
of descent, extreme vibrations, and temporary loss of 
control. Flight profiles where the vortex ring state is 
typically encountered include steep approaches, approaches 
with a tail wind, or conditions where engine power is 
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limited. The vortex ring state can also be encountered 
during aggressive maneuvering flight where flow velocities 
are normal to the rotor disk plane. 

Among different aviation communities the term vortex 
ring state can be referred to by several different 

meanings. In Naval Aviation, the vortex ring state is 
referred to as "power settling." This phenomenon should 

not be confused with the flight profile "settling with 
power," which often refers to situations where power 
required exceeds power available. However, this condition 
can evolve into encounters with the vortex ring state. When 
used in this thesis, the term power settling refers to the 
vortex ring state. 

3. The Windmill Brake State - Autorotation 

As the rate of descent is increased, the flow up 
through the rotor increases. When the value of this flow up 

through the rotor reaches a value twice that of the hover 

induced velocity, the rotor is considered to be in the 
windmill brake state. The rotor is now extracting energy 
from the high velocity slipstream. The rotor disk can be 
thought of as an effective parachute. Figure 1. In fact, 
Gessow and Myers [Gessow, 1952] show that the lift force 
generated by a rotor in autorotation is equivalent to that 
produced by a parachute of the same diameter as the rotor. 
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II. ANALYTICAL AND EXPERIMENTAL MODELS OF THE 

VORTEX RING STATE 


A. FUNDAMENTAL MOMENTUM THEORY 

Momentum theory and blade element theory are used to 
analyze rotor performance in a hover, climb and descent. 
Momentum theory applies the basic conservation laws of 
mass, momentum, and energy. Blade element theory predicts 
the actual forces acting on an element of the rotor blade 
as a result of its motion through the air. Figure 3 
displays the basic steady flow parameters and the funneling 
effect. 



Figure 3. Flow Field Condition. [From: Prouty, 1986] 


7 




The rotor is modeled as an actuator disk with zero 
thickness that supports a pressure differential and 
accelerates air. The velocity is located at the rotor 

plane and Vj is at a point below the rotor, as in Figure 3. 
Mass flow through the rotor system is given as: 

m = pAVj 

Where p is the density of air (associated constant) and A 
is the area swept out by the rotors. The thrust T produced 

by a rotor system is a product of the mass flow m, and the 
change of velocity of the air mass AV . 

T = mAV 

Assuming the air far above the rotor Vq has no velocity, AV 
is simply V 2 . Then rotor thrust is given as: 

T = mV2 = pAVy2 

Kinetic energy, KE, is imparted as the air mass passes 
through the rotor. 

EnergyRate = Force{T)xVelocity(y 1) 

The change in kinetic energy is the rate energy is 
imparted to the downstream wake. The change in kinetic 
energy across the entire system is simply the kinetic 
energy at region 2. The rate at which energy is imparted to 
the flow stream by the rotor must equal the rate at which 
kinetic energy is generated in the wake. 



dKE 

dt 


J wake 


1 

2 


mV^ 


Setting the two equations equal: 




Considering a closed system, the mass flows are equal, 
due to continuity. Leaving the following expression: 

Substituting this back into the thrust equation 
reveals the following equation for hover induced velocity: 


Vn = 



This equation has been derived assuming a uniform 

flow field cross section. Rotational energy induced by 
torque from the rotor is neglected. The purpose of blade 
twist is to produce a more uniform flow field across the 
rotor disk. Anything other than uniform flow through the 
rotor plane increases the kinetic energy for the same 
momentum value and is less efficient. 

1. Vortex Ring State 

Momentum theory remains valid until the flow field 
begins to breakdown and become unsteady, as in VRS. This 

can occur at a moderate rate of descent, 500-700 feet per 

minute. The flow in the recirculated wake becomes difficult 
to analyze. When the rate of descent reaches the hover 
induced velocity there is analytically no mass flow through 
the rotor. Because of this fact, empirical data from 

instrumented rotor systems is used to obtain performance 
results. Momentum theory no longer is valid during 
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operation in the vortex ring state because of the ambiguity 
in defining a control volume where the governing equations 
can be applied. 

2. Autorotation 

Autorotation is considered to occur at the boundary of 
the vortex ring state. That is, ideal autorotation occurs 
when the rate of descent equals twice the hover induced 
velocity. During power-off descent the source of power 
rotating the rotorhead is the air flowing upward. The lift 
vector of the blade element is equal to the drag component. 
Using blade element and momentum theory, an approximation 
for rate of descent during an autorotation is Vo~2[V^]. 

3. Windmill Brake State 

At very large rates of descent the rotor is extracting 
power from the air flowing upwards through the rotor. 
Momentum theory can be applied to determine rotor 
performance when operating in the windmill brake state with 
a clearly formed slipstream. 

B. WOLKOVITCH 

Dr. Julian Wolkovitch published a paper in the Journal 
of the American Helicopter Society in 1972 entitled 
"Analytical Prediction of Vortex Ring State Boundaries for 
Helicopters in Steep Descents." He used momentum theory, 
supported by empirical data, to predict the upper and lower 
boundaries of the vortex ring state. 

The vortex that forms around the circumference of the 
rotor disk during a vertical descent was used to calculate 
an average value for the velocity located at the core of 
the vortex. Wolkovitch defined the upper boundary of the 
vortex ring state when this resulting velocity, acting at 

the core, equaled zero. When the vortex core has no 

10 



movement away from the rotor plane, unsteady flow will 
occur. This was calculated to be the critical descent 
velocity () at which the vortex ring state occurred. 

The critical descent velocity was also greatly 
dependent upon rotor blade tip loss that occurred due to 
non-uniform lift distribution across the blade. Lift 
rapidly decreases across the span near the tip of the blade 
during encounters with the vortex ring state. This results 
in a corresponding loss of induced flow through the rotor 
and decreases the effective radius of the rotor system. 

To substantiate his results, Wolkovitch used 
experimental data showing thrust fluctuations indicating 
unsteady flow when operating in the vortex ring state. His 
results showed significant thrust fluctuations at 
approximately Krir=0-7V'^- He reported mean-to-peak thrust 
fluctuations of 14 percent of the mean thrust, as shown in 
Figure 4. 
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Figure 4. Measured Thrust Fluctuations During a 
Vertical Descent. [From: Wolkovitch, 1972] 

The windmill brake state occurs when the descent rate 
is increased to a point where the steady flow reappears. 
This defines the lower boundary of the vortex ring state. 
The entire rotor wake now exists above the rotor plane. 
Wolkovitch suggest that the lower boundary is not as well 
defined as the upper and therefore uses an empirical factor 
to account for the contraction of the wake above the rotor 
during descent. 

Wolkovitch compares his calculated boundaries to that 
obtained from flow visualization. Figure 5. It is apparent 
that the upper and lower boundaries do not merge as 
airspeed is increased, implying that vortex ring state can 
be encountered at high airspeeds. 
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Figure 5. Calculated VRS Boundaries Compared to Flow 
Visualization Boundaries. [From: Wolkovitch, 

1972] 


C. PETERS AND CHEN 

Drs. David Peters and Shyi-Yaung Chen published a 
paper in the Journal of the American Helicopter Society in 
1982 following up on Wolkovitch's work but using a more 
consistent wake model. [Peters, 1982] They refer to the 
hover condition when the rotor's induced flow is imparting 
energy in to the flow field and the windmill brake state as 
extracting energy. 

They took into consideration the horizontal component 
of velocity generated by forward flight. They derived an 
equation that relates the normalized rate of descent T] , 

horizontal velocity |I , and hover induced flow . 
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Peters and Chen conclude that during high values of 
inplane flow, the vortex ring state does not occur. 
[Peters, 1982] They conclude that there is a vortex ring 
state boundary at a determined forward airspeed specified 
by a maximum normalized inplane horizontal velocity. 

From actual flight data there is clearly a maximum 
inplane velocity for which the vortex ring state 
dissipates. Increasing forward airspeed is a typical 
recovery technique from the vortex ring state. The unsteady 
air formed around the rotor tips is swept aft and above the 
aircraft. With consideration of inplane velocity, the 
subsequent equations give rise to a more defined boundary 
for the vortex ring state. Figure 6. 
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Inplane Velocity 

0.2 0.4 0.6 0.8 



Contours Of Vortex Strength 
Peters and Chen Boundary 

Figure 6. Peters and Chen VRS Boundary. [From: 

Peters, 1982] 

Peters and Chen showed that there is a boundary for 
the vortex ring state at a normalized inplane forward 
velocity and proved the general boundary was larger than 
predicted by Wolkovitch. They also determined that the 
center, not the edge of the region is located at normalized 
rate of descent of 0.707. 

D. GAO AND XIN 

Gao and Xin published a paper in the First Russian 
Helicopter Society Annual Forum Proceedings in 1994 
entitled "An Experimental Investigation into the Vortex 
Ring State Boundary." 
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Goa and Xin's primary concern was to eliminate the 
common wind tunnel limitations, interference flow 
interaction at the tunnel walls causing boundary layer 
interference, and the poor low speed characterists of wind 
tunnels. The disturbed air surrounding the rotor during the 
vortex ring state can extend out a distance of several 
rotor diameters. Therefore, the interaction with the tunnel 
walls can be substantial. 

They used a whirling beam device capable of measuring 
torque and thrust with a strain gauge. Output from the 
strain gauge was amplified and sampled. They tested three 
different airfoils with varying twist and chord. 

Gao and Xin used a semiemperical approach and 
momentum theory to determine new boundaries of the vortex 
ring state. They experimentally determined that the 
boundary occurs when the free stream velocity reaches a 
critical value of 0.28V';,. Descent velocities 0.6 to 0.8 
times the hover induced velocity ( 0.6 V;, < V < 0.8 V;, ) is 
identified to be the region of most sever encounters with 
VRS. Data indicated drastic fluctuations in torque and 
thrust during operation in the VRS region. Figures 7 and 8. 
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Figure 7. Torque Fluctuations During Vertical 
Descent. [From: Gao And Xin, 1994] 

Amplitude Of Thrust Fluctuation/Mean Thrust % 
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Figure 8. Thrust Fluctuations During Vertical 
Descent. [From: Gao And Xin, 1994] 
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A significant discovery was the correlation of mean 
torque to descent rate; they called this the power settling 
phenomenon. The "bucket" torque setting can easily be seen 
in Figure 9 at a descent rate of V=0.28 . Beyond this 
point it requires an increase in mean torque to descend at 
a faster rate. Above Gao and Xin predict the 
increase is due to local blade stall. Although, above 
descent rates greater than V=0.the torque fluctuations 
decrease as the rotor system transitions into the windmill 
brake state. They found that the lowest descent rate where 
the vortex ring state can be encounter is V=0.2. 
Therefore, this value is defined as the critical descent 
velocity, Vcrit=0.28 . 


Mean Torque/Mean Torque In Hover 



Figure 9. Mean Torque During a Vertical Descent. 

[From: Gao And Xin, 1994] 
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Gao and Xin determined for their rotor system, that 
the vortex ring state is encountered only for steep 
approaches of greater than 70 degrees. 

All three models are displayed in Figure 10. 



♦ Gao and Xin Boundary 
■ Peters And Chen Boundary 
A Wolkovitch Lower Boundary 
Wolkovtich Upper Boundary 
X Flow Visualization 

Poly. (Wolkovtich Upper Boundary) 
Poly. (Wolkovitch Lower Boundary) 


Figure 10. Superposition of Vortex Ring State 

Boundaries and Flow Visualization. [From: 
Varnes, 1999] 


E. ENGINEERING MODEL 

Several recent analytical models for predicting the 
flight boundary of the vortex ring state have been 
developed. An engineering model was developed based on 
defining VRS development in terms of the rate of growth of 
the disturbances relative to the rate at which disturbances 

are covected downstream [Newman, 2001] . The growth of 
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disturbances in the flow field is proportional to the wake 
vorticity. The shed vortices are a direct result of rotor 
thrust. The engineering model uses conventional actuator 
disk theory to identify a rate of transport of the wake 
instabilities into the downstream flow field. The 
foundation of the theory is based on classical momentum 
theory. [Newman, 2001] 

F. WAKE STABILITY THEORY 

Using a linear analysis of the rate of growth of 
perturbations introduced into the wake structure can help 
identify the process of degeneration. The wake stability 
theory assumes an incompressible, irrotational fluid with 
embedded vorticity which will be transferred to the 
surrounding fluid. The stability of the rotor wake is 
dependent on both rotor geometry and the operating state. 
With this method, predicting what happens after the rotor 
wake has become unstable is not possible. Numerical methods 
must be used to reveal the non-linear formation of the 
unstable wake. [Brown, 2002] 

G. VORTICITY TRANSPORT MODEL (VTM) 

The vorticity Transport Model (VTM) developed by Brown 
uses a direct numerical computational solution to predict 
the development of the rotor wake vorticity. Using the 
Navier-Stokes Equation, assuming incompressibility and the 
limit of vanishing viscosity, shows rotor wake vorticity is 
related directly to aerodynamic loads on the blades. The 
Biot-Savart relationship is used to relate the velocity at 
any point near the rotor to the vorticity distribution. 
Aerodynamic loading on the airfoil is determined from the 
geometry and strength of the rotor wake. Moderate thrust 
fluctuations at low descent rates quickly increase as the 
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descent rate approaches 0.5 times the hover induced flow. 
[Brown, 2000] The model can represent blade-to-wake, and 
wake-to-wake interactions that lead to growth and rupture 
of the vertical wake pattern. [Brown, 2002] 

H. FREE VORTEX MODEL (FVM) 

The FVM is a computational method with improved 
computational efficiency, compared to VTM. Although, there 


is a 

loss 

of 

resolution when using 

FVM. 

A 

Lagranian 

description 

of 

flow is 

used 

to pattern 

the 

FVM. 

Discrete 

line 

vortices 

represent structure of 

the 

rotor wake 

vorticity. 

For 

straight 

line 

segmentation 

of 

the 

vorticity 

being 

shed 

and 

trailed 

from 

the blades 

the 

approach is 


spatially second order accurate. [Brown, 2002] 
Unfortunately, computational model descriptions of highly 
unstable, nonlinear, time dependent flow in the vortex ring 
state are susceptible to small perturbations and can result 
in rapid divergence of FVM and VTM predictions. 
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Ill. EXPERIMENTAL APPARATUS 


The analytical and experimental research presented has 
provided extensive understanding and further defined the 
boundaries of the VRS. Although published research of 
studies conducted in a water tunnel is very limited. It is 
believed that vibrational analysis of a rotor operating in 
VRS deserves additional investigation. In the effort of 
contributing to this field of research a small scale rotor 
was constructed to operate in a water tunnel. The apparent 
lack of vibrational analysis of the rotor operating in VRS 
has motivated this thesis. 

A. WATER TUNNEL 

A flow visualization water tunnel manufactured by 
Eidetics International was used to conduct our experiment. 
Figure 11. The tunnel has a 15 inch by 15 inch flow area 
and can produce smooth flow at speeds up to 16 inches per 
second. 
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Figure 11. Flow Visualization Water Tunnel with 

Scaled Rotor Installed. 


The advantage in using a water tunnel for this 
research is the ability to reduce the flow velocities 
compared to full scale models operating in air. The 
similarity of flow dynamics between water and air enable 
this advantage; the streamline patterns are geometrically 
similar and the force coefficients are the same. 
Distributions of parameters such as VIV^ are the same when 
plotted against common nondimensional coordinates. The flow 
over geometrically similar bodies at the same Mach and 
Reynolds numbers are dynamically similar, therefore the 
lift, drag, and moment coefficients are identical. 
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Disregarding the scaling of rotor chord length, c, 
comparison of Reynolds numbers associated with water and 
air reveals the source of the flow velocity advantage. 
Reynolds number calculations indicate a reduction of 
velocity by a factor of 15. p is density and |I is the 
dynamic viscosity. 


R. = 


pVc 

P- 


(^) =(hi) 

''water ^ ''air 

|i 

V = (-Hi) . y. 

P |i 

V =—V ■ 

water ^ ^ air 


Rotor tip speeds of a full scale helicopter are 
approximately 600 ft/sec. The chord length can range from 
six to eighteen inches. Therefore, Reynolds numbers of a 
full scale rotor blade operating in air are on the order of 
10*^. 

The chord length of the scaled model blade used in the 
present experiment is one inch. A blade tip speed of 14 
ft/sec was chosen. This results in a Reynolds number in 
water on the order of 10^, close to that of a full scale 
blade. In addition, the corresponding rotational speed 
provided excellent flow visualization. 

B. SCALED ROTOR SYSTEM MODEL 

A scaled model rotor system was constructed using the 
tail rotor system of a Hawk IV radio controlled hobby 
helicopter. Figure 12. The rotor diameter is nine inches 
and corresponding rotor area is 0.44 ft^. A variable RPM 
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Bodine electric motor provided 1/8 horsepower to rotate the 
two symmetric blades at the desired speed of 360 RPM. The 

blade pitch angle was set at ten degrees of pitch and 

tunnel flow velocity was varied from one to nine in per 
second. 

Small compressed air bubbles were injected into the 
flow field at the edges of the rotor plane. The air 

discharge nozzle was crimped to produce as small bubbles as 
possible. The viscosity force acting to hold the bubbles in 
solution countered the buoyancy force and helped keep the 
air bubbles from being forced to the surface. Two air 

outlet probes were placed at points with very close 
clearance to the rotating blades. Each passing blade picked 
up the injected air and formed a flow field of small 
bubbles. 



Tank Flow 
Direction 


Load Ceil, 
Wheatsone Bridge 


Figure 12. 


Rotor System Device. 
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C. DATA ACQUISITION 

For static measurements a Wheatstone bridge strain 
gauge load cell was connected to a meter, and provide a 
visual output display of the thrust load. The load cell was 
attached as shown in Figure 12. Calibration was 

accomplished with static weights prior to each experiment. 

The same Wheatstone bridge strain gauge was used to 

measure dynamic forces on the rotor system. The load cell 
signal was feed to an Ectron model 563H amplifier. 5 Vdc 
excitation voltage was provided by the amplifier to the 
Wheatstone bridge. A Hewlett Packard 35665A Dynamic Signal 
Analyzer was then used to compute power spectra from the 
thrust data. Figure 13. All data were low pass filtered at 
bandwidth 64 Hz and sampled at 128 samples per second for 2 
seconds. Four time records were (rms) averaged with 0% 

overlap. The display was 0 to 50 Hz with a frequency 

resolution of 0.5 Hz. A Hanning window was used. The 
Nyquist criteria was met by sampling more than twice as 
fast as the highest frequency component displayed on the 
analyzer. 
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Figure 13. 



Ectron 563H Amplifier and Hewlett Packard 
35665A Dynamic Signal Analyzer. 
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IV. EXPERIMENTAL RESULTS 


A key factor in accurately predicting the aerodynamic 
behavior of a rotor system is an improved understanding of 
the wake distribution generated by the rotor system. The 
understanding and eventual prediction of wake dynamics 
during descending flight is key to improved understanding 
of the phenomenon of the vortex ring state. 

A. HOVERING CONDITIONS 

A helical vortex wake can be seen forming downstream 
of the rotor blades in Figure 14. The rotor system is 
producing thrust to the right and can be visualized as a 
helicopter in flight by rotating the image 90 degrees 
counterclockwise. Images were captured using a digital 
still camera with a fast shutter speed setting of 1/320 
sec. When the flow in the tunnel is turned off, the flow 
visualization section becomes a closed system with no area 
for rotor wash to dissipate. This made simulating hover 
conditions impossible without a slight flow velocity 
established in the tunnel. When a slight tunnel flow 
velocity of three inches per second was turned on, the 
rotor wash was able to displace naturally. The rotor device 
was installed during this test so the flow was from above, 
if the rotor device was viewed as a hovering aircraft (from 
the right of the rotor in the photos) . This can be thought 
of as a climb profile. 
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Figure 14. Helical Vortex Wake in a Climb. Note Wake 

Contraction. 

1. Momentum Theory Validates Observation 

In Figure 14 the rotor was operated at 360 RPM (6 Hz) 
with blade pitch set at ten degrees. In addition to the 
rotor producing thrust, and shedding the helical wake 
downstream, an additional water tunnel flow rate exist from 
right to left in the photo. The velocity, , in the wake 

at a position approximately one inch downstream of the 
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rotor was obtained from the photograph by 

Ar 2in , 

V i. = — =-= 12 m/sec . 

A/ 0.17 r where Ar is the measured distance 

from the bubbles entrained in the core of a helical vortex 
at a position corresponding to one cycle elapsed and A/ is 
the time for one cycle of rotation. This velocity equals 
the sum of the climb induced velocity through the rotor, 
11.^, and the water tunnel flow velocity, 11 =3 m/sec. 

The tunnel flow enabled the rotor wash to be convected 
downstream, and to establish a well-defined flow field such 
as shown in Figure 1, and for which momentum theory 
applies. The tunnel flow velocity of 3 in/sec can be 
thought of as a rate of climb, 11 . 

A load cell thrust measurement {T ) during this 
observation was recorded as 0.5 lb. A hover induced 
velocity (I 4 ) is calculated for this specific thrust value: 



= 6.5 in I sec 

Using momentum theory the following equation gives the 
induced velocity for this climb condition. [Prouty, 1986] 
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= 5.0 in I sec 


The velocity through the rotor is, then +V'^ =8/n/sec . 

Conservation of mass between the rotor and the wake 

downstream requires • Taking = 0.44/?" 

2 

and (Vi^+V^) = S inlsec = — Jils the left hand side is calculated 

to be approximately 0.29 ft^/sec. 

The diameter of the wake where the velocity was 

estimated is estimated to be approximately 7.5 - 8.0 

inches. This results in an estimated wake area = 0.33/?" . 

With the estimated value of V„^i^ = \lin!see, The right hand 
side of the mass conservation equation is estimated to be 
approximately 0.33 ft^/sec, a difference of only 10% from 

the left hand side. This should be considered a remarkable 
agreement, and shows that momentum theory can be used to 
adequately describe the gross features of the flow in this 
experiment. 

Momentum theory predicts a rapid wake contraction 
downstream of the rotor, which is clearly apparent in 
Figure 14. Theory further predicts that the area of the 
wake very far downstream of the rotor should be the rotor 
swept area. From measurements taken off the photograph of 
Figure 14, we observe that, at a position two inches 
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downstream of the rotor, the diameter of the wake is 
approximately seven inches, and the wake area relative to 

(if 

the rotor swept area is already -^ = 0.6 

(9)" 

B. VORTEX RING STATE 

Throughout the experiment the rotor was operated at a 
constant 360 RPM with a blade pitch angle of ten degrees. 
With these parameters, the hover induced velocity was 
estimated by visually observing when the helical vortex 
wakes collapse in the plane of the rotor. The water tunnel 
flow velocity was increased until this occurred. Figure 15. 



•S' . 


Figure 15. Collapsing Vortices, Impending VRS with 

Flow Rate at 4.5 In/Sec. 
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The average thrust measurement when the vortex wake 
collapsed was recorded as 0.7 lbs. A corresponding velocity 
is calculated which can be considered as a hover induced 
velocity with a thrust measurement of 0.7 lbs. 


K=. 


'2pA 


OJlb 


2(1.95"^)7r(0.375/0" 


= 4 . 5 ^ 
sec 


The collapsing of the helical vortex wake structure 
into the plane of the rotor indicates the onset of the 
vortex ring state. Recirulation of the rotor wake began at 
this point. Figure 16. The test also showed some water 
tunnel "wall effects" due to constraints of the small flow 
section of the water tunnel. 



Figure 16. VRS Recirulation Take at Slow Shutter 

Speed. 
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C. TRANSITION TO AUTOROTATION 

As the flow velocity of the tunnel was increased 
beyond 4.5 in/sec the vortices reformed on the opposite 
side of the rotor, the upstream side. There was a clear 
transition from VRS into autorotation. A helical vortex 
wake flow pattern, similar to that seen in a hover, 
reformed on the opposite side (upstream) of the rotor 
system. This would compare to the wake forming above a 
helicopter during an autorotation. A clear transition from 
hover through VRS and then into autorotation was achieved 
by gradually increasing the water flow velocity in the 
tunnel. 

D. VIBRATION ANALYSIS 

In flight, the primary indicator of the vortex ring 
state in flight is the associated increase in airframe 
vibration levels. Erratic thrust fluctuations and irregular 
blade flapping cause the increased airframe vibrations. 
Once fully operating in the vortex ring state the vibration 
levels are the most severe encountered in any flight 

regime. 

A paper published by John E. Yeates in 1958 at NASA 
Langley explored the vibrational effects and subsequent 
rotor blade stresses produced during operation in the 

vortex ring state. A tandem rotor aircraft, with three 

blades on each rotor, was used in the testing. Sensors 
attached on the front and rear masts measured vertical 
vibrations. 95% of the vibratory response was reported at 
three times the operating frequency (3P) . Although, 

instrumentation was not onboard to record individual 
harmonics during operation in the vortex ring state. [Ref. 
Yeates, 1958] 
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The vortex ring state can be viewed as a form of wake 
instability. In a descent, the random aerodynamic 
excitation builds up. As the instabilities reach the rotor 
plane they begin to effect the aerodynamic section loading 
across each airfoil. Operation in this adverse condition of 
VRS does not allow the downwash to convect away from the 
rotor. Instead, the vortex wake accumulates in the plane of 
the rotor. As blade vortex interaction repeatedly occurs 
due to unsteady blade air loads, the result is increased 
airframe vibrations. 

1. Experimental Vibration Results 

It was found that the increase in amplitude of higher 
harmonic frequency variations of thrust were a good 
indicator of impending vortex ring state. In this 
experiment the rotor operating speed was 6 Hz (IP), on a 
two-bladed rotor system. As the vortex ring state was 
encountered, higher harmonics of the blade passage 
frequency began to increase in amplitude. The content of 
the thrust frequency input at higher than n/rev (where n 
equals the number of rotor blades) is found to be a 
significant feature to the vibration spectrum associated 
with operation in the vortex ring state. "One per rev" (IP) 
content of thrust variation is primarily driven by cyclic 
pitch and resultant blade flapping. It was found to be 
relatively invariant to impending VRS and therefore is not 
a good indicator. 

In preparation for data collection, a vibration shaker 
test was conducted on the rotor device to determine its 
natural frequencies. Results of the shaker test indicated 
that natural resonant frequencies appeared to be above the 
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test bandwidth of 50 Hz. Although this test was not 
conducted in the water, it is believed that conditions when 
in the water would tend to dampen natural resonances. 

The experiment was conducted while operating the rotor 
at a constant 6 Hz rotational rate. The flow in the tunnel 
was gradually increased to 4.5 inches per second, where the 
helical vortex wake collapsed into the plane of the rotor. 
Continual operation in this state was possible. The 
visually recirculating flow field was a clear indication of 
operation in the vortex ring state. 

Thrust data recorded during VRS, presented in the 
frequency domain, clearly indicates the increase in 
strength of higher harmonics. The blade passage frequency 
is defined as n times the operating frequency, P, where n 
is the number of blades, i.e., 2P (12 Hz) in the case of a 
two bladed rotor. The even multiples of the blade passage 
frequency, 4P (24 Hz) and 6P (36 Hz), are shown to increase 
during VRS, Figure 17. This response is expected as blade 
vortex interaction increases, causing increased variations 
in the section aerodynamic loading. 
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Figure 17. Scaled Vibration Spectrum from Hover (0 
In/Sec), VRS (4.5 In/Sec) and Autorotation 
(9 In/Sec). All Are Plotted Against Same 
Vertical Scale, Which Is Normalized to the 
6 Hz Autorotation Peak. 
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Note that significant hub vibration response only 
occurs for frequencies in the vicinity of the rotor 
rotation frequency (IP) and at even multiples of the blade 
passage frequency (2P). This is expected behavior. Rotating 
blades produce a vertical force at the root of the blade 
which results in a vertical excitation of the hub. IP (one 
times rotation frequency) is primarily driven by cyclic 
pitch or can be a result of out of balance blades, most 
likely the cause in our data. Assuming all blades see the 
same loading at the same azimuth, only even multiples of nP 
frequency will be transmitted through to the fixed system, 
the fuselage, or hub attached to the support structure in 
our case. [Prouty, 1986] 

Vibrations of the airframe at even multiples of n/rev 
are a result of (n-1), n and (n+l)/rev loads in the rotor 
system, which are filtered through the rotor as n/rev. (n) 
is the number of blades on the rotor. Harmonics of the 
airfoil loading cause this response. This corresponds to 
2P, 4P, and 6P present in the fuselage of a two blade 
rotor. For example, the 4P frequency in the fuselage, or 
hub, is a result of 3P, 4P and 5P loads in the rotor 
system. 

The splitting in the peaks in the hub vibration power 
spectra (Figure 17) is not expected, and the cause is not 
understood. It appears to be of the form of an amplitude 

modulation cos((0^gjO cos((00 = ^[cos((0+ (0^^,^)?+cos((0-(0^^,^)?] . This is 

indicative of a primary frequency, CO , operating in 
conjunction with an additional modulation frequency, (Oj„od/■ 
superimposed. Therefore, in quantitatively reporting the 
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results of the power spectrum analysis. Table 1, the 
amplitudes of the split peaks were combined as the square 
root of the sum of the squares. Peak values at 2P, 4P, and 

6P are reported, and the higher harmonics (4P, 6P) are 
compared to values at 2P. 



Hover 

VRS 

Autorotation 

2P (averaged) 

0.36 

0.38 

0.35 

4P (averaged) 

0.26 

0.33 

0.14 

6P (averaged) 

0.28 

0.46 

0.17 

4P/2P 

0.72 

0.87 

0.40 

6P/2P 

0.78 

1.21 

0.49 


Table 1. Maximum Amplitudes and Comparison of Higher 

Harmonic Magnitudes of Harmonics Listed in Figure 17. 

The data collected during hover condition are believed 
to be affected by the confines of the water tunnel. The 
inadequate dissipation of the rotor wash may have created 
greater blade vorticity interaction and some increase in 
higher harmonic vibration levels than would free field 
conditions. Even so, considering this result, the 4P and 6P 
harmonic are clearly higher during VRS. 

The flow velocity was increased to nine inches per 
second to achieve windmill brake flow conditions. The 
helical vortex wake flow pattern reformed on the opposite 
side of the rotor system. Clear smooth flow patterns were 


40 




observed at this flow velocity of approximately two times 
hover induced, 2V'^ . The higher harmonics at 4P and 6P 
decreased in amplitude. 

It should be noted that the flight profile of a full 
scale helicopter does not correlate exactly with this 
experiment. Autorotation in a full scale helicopter is 
initiated by sharply reducing the blade pitch angle 
(collective). In our experiment the operating RPM and blade 
pitch remain constant. As the tunnel flow velocity is 
increased, this increases the angle of attack and thrust 
loading. Data recorded during an actual autorotation would 
be associated with a much lower thrust loading than our 
corresponding data. Therefore, the amount of higher 
harmonics observed in the present experiment under 
autorotation conditions is probably relatively greater than 
would be observed under full scale flight conditions, 
compared to the VRS state. 

Corresponding increase in higher harmonics during VRS 
would further increase if operated under full scale flight 
condition. 

To correlate these model test results, an 
investigation was carried out with respect to full scale 
data collected by NASA in 1964. [Scheiman, 1964] Under the 
direction of James Scheiman, an extensive amount of in 
flight data were collected on a fully instrumented H-34 
helicopter. These are the only full scale data ever taken 
where a significant number of VRS conditions were recorded. 
An analysis of the H-34 NASA data revealed similar results 
as recorded in the water tunnel experiment and are 
discussed in a later section. 


41 



2. Karman Vortex Street Effect 

The Karman vortex street effect is here in considered, 
to determine its potential effect on the vibration data 
collected. The Karman vortex street is a regular array of 
vortices of alternating vorticity, convected downstream of 
a blunt obstacle. Figure 18. A Karman vortex street formed 
around the long drive shaft and transmitted transverse 
vibrations into the apparatus. Although the vibration 
amplitude from the Karman vortex street was assumed to be 
minimal, an analysis was conducted. 



Figure 18. Karman Vortex Street Around a Circular 
Cylinder. [From: Van Dyke, 1982] 


The Reynolds number associated with the drive shaft at 
a tank flow speed {V ) of 4.5 in/sec is calculated, (c) is 
the diameter of the cowling around the drive shaft. 


^ (0.036^)(4.5—)(0.75m) 

P _ in sec 


R 




55.89;cl0 


.-6 


lb 


2174 


m - sec 


With = 2174 , turbulence develops soon after boundary 
layer separation, but the boundary layer is still laminar 
around the cylindrical shaft. Strouhal number is used to 
calculate the frequency (/) of turbulence departing the 
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vertical drive shaft, for a Reynolds number of 2174 the 
value of the Strouhal number is about 0.2. The diameter of 
the shaft {d) is 0.75 inches and {V ) is the free stream 
flow velocity. 


/ 


5 ,= — - 0.2 

' V 

(4.5 m/sec)0.2 
OJ Sin 


For 10^ < Re < lO'^ 


1.2 Hz 


The frequency of induced vibrations due to formation 
and detachment of vortices is expected to be approximately 
1.2 Hz. 


Transverse shaft oscillations induced by the regular 
shedding of vortices was considered as a potential cause of 
amplitude modulation of the thrust variations. If this were 
the cause, however, we should expect to observe an 
increasing frequency splitting with increasing tunnel flow 
velocity, which we do not. Hence, we conclude that the 
effects of vortex shedding on the vibration spectrum are 
not significant. 

E. AIR BODY FORMATION UNDER ROTOR DURING VRS 

A peculiar phenomenon occurs as a helicopter rotor 
system continuously circulates air when experiencing the 
vortex ring state. A bubble of air forms under the rotor, 
which appears to grow and burst on a periodic basis, about 
once every second. [Prouty, 1986] This periodic dissipation 
causes a large scale disturbance in the surrounding flow 
field. The bubble erupts, causing erratic flapping in pitch 
and roll as rotor blades interact with this disturbed air. 
This unsteady turbulent air can persist until the rate of 
descent reaches values approximately twice the hover 
induced velocity. 
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The classic Drees and Hendal [Drees, 1951] smoke video 
from 1951 shows an air body that appears to form under the 
rotor during operation in the vortex ring state. Figure 19. 
This air body seems to accumulate energy and then 
periodically burst at intervals of one to two seconds. 
Reports in the literature on VRS comment that the bursting 
interval of approximately one second appears to be 
invariant to rotor size. That is, the bursting period is 
about the same for rotors regardless of size from model 
rotors up to full scale rotors. 



Figure 19. Vortex Ring State Air Body Formation. 

[From: Prouty, 1986] 
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The power spectrum data recorded during Goa and Xin's 
experiment indicated an aerodynamic fluctuation on the 
order of 1 Hz, Figure 20. This appears to coincide with the 
bursting of the air body set up under the rotor during 

operation in the VRS. 

70 
60 
60 
40 
CJO 
20 
10 
0 

0 2 4 6 8 10 ta !4 16 18 20 22 24 

Frequency (Hz) 

Figure 20. Power Spectrum During a Descent. [From: 

Gao and Xin, 1994] 

During VRS conditions in the water tunnel, the same 

flow manifestation was observed. The phenomena manifested 
itself by periodic large eruptions occurring at the surface 
of the water at a periodic rate of 1-2 Hz, Figure 21. The 
formation of this air body, and its periodic dissipation, 
appear as though it may be a three dimensional Karman 
vortex street effect caused by the rotor disk acting as a 

solid body in the flow stream. 
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Figure 21. Periodic Surface Disturbance During VRS. 

A contrasting view of normal smooth flow is during the 
autorotation profile as shown in Figure 22 . It is clear 
that the periodic disturbance at the surface is no longer 
present. 

Note, the small dip on the top right surface of the 
water in Figure 22 is due to the incompressibility of water 
and the tight operating constraints of tunnel flow area. As 
flow is accelerated through the rotor system, the surface 
is sucked down. During autorotation the rate of descent 
equals twice the hover induced velocity. The flow pattern 
is convected around the rotor and around the disk, causing 
a corresponding dip on the surface. 


46 






Figure 22. Smooth Flow Field. 


F. V-22 BLADES 

NAVAIR researcher Kurt Long at NASA Ames generously 
provided blades with twist similar to those on the V-22 . 
The symmetric blades used up to this point had zero twist 
(i.e., constant pitch) . The V-22 blades have a high degree 
of twist at the root of the blade and transition to zero 
twist at the tip. These blades were installed and tested 
under the same parameters as used for the symmetric blades. 
Thrust variation along the rotor blade, due to its pitch 
variation, seems to affect the hover induced velocity, as 
observed by the spacing of the helical vortex wake. The V- 
22 blades, with high twist, showed counterintuitive 
differences when subjected to operation in the vortex ring 
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state. Wake stability, affected by the thrust distribution 
along the blade appears to have a large effect on the 
descent rate boundary at which VRS is encountered. With 
higher rotor disk loading the boundary occurs at higher 
rates of descent than for a symmetric airfoil. Extensive 
testing of the V-22 blades was limited due to time 
constraints of this project but further investigation using 
the water tunnel promises to produce interesting new 
results. 
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V. 


SCHEIMAN H-34 FLIGHT TEST DATA 


A. FLIGHT DATA ANALYZED 

Under the direction of James Scheiman, in 1964, NASA 
Langley Research Center instrumented an Army Sikorsky H-34 
helicopter to obtain rotor blade load distribution, 
flapwise bending moments, chordwise bending moments, 
harmonic analysis of blade root motion, and section 
aerodynamic loading. [Scheiman, 1964] Flight data are given 
for hover, climbs, autorotation, and of particular 
interest, partial rate of descent at slow airspeed. 
Additionally, there are insightful corresponding pilot 
comments. A comparison was made of the harmonic analysis 
presented in Scheiman's report with the harmonic analysis 
of the present experiment. 

In Scheiman's study the aircraft had four rotor blades 
and operated at the following conditions during the 
selected flights listed in Table 2. Approximate gross 
weight was 11,500 lbs and air density was approximately 
0.00216 slugs/ft^. 
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Forward 


Vd, av 


Pilot Remarks 


Table Flight 


Power 


Airspeed ^ 


(in Hg) (ft/min) 


75 

56 

0 

31.4 

1200 

Rough; blades flapping 

erratically; unsteady flight 

79 

59 

8 

39.3 

2100 

Helicopter temporarily out of 

control 

80 

59 

7 

40.1 

1850 

Helicopter temporarily out of 

control 

84 

62 

6 

37.4 

2600 

Unsteady flight 

4 

01 

0 

38.5 

0 

Hovering in light wind 


Table 2. Scheiman - NASA Test Data. 


The first four flights listed in Table 2 were 
specifically chosen because of the clear indications of 
operation near the vortex ring state. The remarks by test 
pilots indicating uncontrolled high descent rates, 
vibrations, and temporary loss of control are all typical 
results of operating in the vortex ring state. These 
comments, in conjunction with the supporting flight test 
data, are clear indications of operation in the vortex ring 
state. 

1. Section Aerodynamic Loading 

Aerodynamic loading data from hover flight are 
compared to data taken during flight in the vortex ring 
state. Flight number 79, encountering a descent rate of 
2100 feet per minute and corresponding comments of 
"helicopter out of control," is clearly encountering VRS. 
Data from flight 4 are from a hover, and provide a baseline 
for comparison. The term r/R is the radial distance along 
the span of the rotor blade, moving radially outward. A 
point at r/R=0.95 is chosen to display the load 
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distribution corresponding to blade azimuth. Figure 23. 
Azimuth, \|/ , is in degrees of blade rotation, referenced to 
zero at a line extending aft along the centerline of the 
tail boom. In a hover, the distribution is smooth, except 
for a dip as the blade passes over the tail boom. In 
contrast, the loading associated with the VRS flight 
conditions exhibits drastic variations, due to blade vortex 
interaction. 



Figure 23. Section Aerodynamic Loading at R/R=0.95 as 
Measured in Hover and VRS at 2100 Fpm Rate 

of Descent. 


2. Flapwise Bending Moment 

Flapwise bending moment data from flights in the VRS 
were compared to data collected during hover. Point r/R = 
0.575 was chosen to reduce some of the higher harmonic 
content resulting from the erratic flapping further out on 
the blade. A comparison of hover data to that from a flight 
in known vortex ring state conditions is shown in Figure 
24. [Varnes, 1999] 
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Figure 24. H-34 Flapwise Bending Moment At R/R=0.575. 


Analysis of the flapwise bending moment data indicate 
highly erratic behavior during VRS conditions. Broadly 
varying flapwise bending moments are evident. Compression 
of the upper surface of the blade is indicated as positive. 
Amplitude changes, as well as sharp variations between 
compression and tension of the upper surface of the blade 
are shown for the case of VRS. Corresponding large blade 
deflection is associated with such values. The drastic 
reversal in compression and tension indicate clearly the 
many blade vortex encounters expected during VRS. Large 
discontinuities are produced, since blade and bending 
moment is the second derivative of blade deflection, the 
rate of change of the blade slope, or gradient. 

3. Harmonic Analysis 

Harmonic analysis of section aerodynamic loading was 
compiled on the five flights listed in Table 2. Considering 
flapwise and chordwise bending, the following equation was 
used to compute section aerodynamic loading. 
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10 

4 (¥) = A) Y, [A cos n(\\f ) + sin n(^^ )] 

n—l 

L^(lblin) 

M^ilblin) 

The section aerodynamic loading, 4(¥) is the integral 
of chordwise differential pressure per unit span (Ib/in) 
over the blade section length. The integrand was 
approximated as a sum over ten stations at equal intervals 
over the full length of the rotor blade. and are 

harmonic series coefficients of section aerodynamic loading 
per unit span, Ib/in. Resulting harmonic amplitudes are 
obtained by calculating the square root of the sum of the 
squares of the coefficients. Figure 25 shows the resulting 
harmonic amplitudes of section aerodynamic loading, and 
compares four flights encountering VRS to data taken in a 
hover. 

In the vortex ring state region extremely abrupt 
changes in airload are a result of blade vortex 
interaction. Higher harmonic content is prevalent during 
this period and is the principal contributor to the 
aircraft vibrations. Study of Figure 25 shows that harmonic 
analysis of the aerodynamic blade loading may be used to 
determine the onset of the vortex ring state. There is 
clearly an indication of the increase in amplitude of 
higher harmonics as the rotor system encounters VRS. The H- 
34 has four blades, and we note the corresponding maximum 
vibration amplitude is seen at 4P, which is four times the 
operating rotor speed, the blade passage frequency. 
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Figure 25. Harmonic Analysis (nP). 


Note that Figure 25 shows significant response at all 
multiples of the rotor operating frequency (IP), whereas 
Figure 17 does not. This is because Figure 25 shows the 
section loading spectrum for a single blade, where Figure 
17 shows the thrust spectrum for the hub, which filters the 
single-blade response. Rotor blades are attached at the 
root to the rotorhead and transmit vibratory responses to 
the fuselage, or hub, as vibratory shears and moments. Flap 
and chordwise bending moments produce vibratory shears at 
the blade root. [Wood, 2000] The rotor head acts like a 
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filter and only allows frequencies of n/rev (the blade 
passage frequency) and even multiples of n/rev to pass 
through. 

The amplitude of the aerodynamic loading gives rise to 
higher harmonics but only 4P and 8P will be transmitted 
through the rotor system into the airframe. Even multiples 
of 4P (i.e., 8P, 16P) will begin to appear in the fuselage 

with increasing amplitude as the onset of increased blade 
vortex interaction occurs. Lower harmonics of the blade 

loading are deemed the most prevalent but higher harmonics 
can be detected. 

A comparison of 4P and 8P harmonics to the blade 

passage frequency (4P) gives a quantifying comparison of 

the amplitude increases. Table 3. These ratio comparisons 
can be used to correlate data recorded during the water 
tunnel investigation. Although, the rotor systems have 
differing numbers of blades, the ratio of the harmonics 

show correlation. Compared to hover conditions, there is a 
significant increase during VRS of the higher harmonic, 8P, 
when compared to the consistent blade passage frequency. 
Although, this does not take into account the increase in 
4P hub vibrations due to 3P and 5P or the corresponding 
increase in 8P due to 7P and 9P. 



Hover 

VRS 

4P (averaged) 

1.8 

2.1 

8P (averaged) 

0.2 

0.7 

8P/4P 

0.1 

0.3 


Table 3. Maximum Amplitudes and Comparison of Higher 

Harmonic Magnitudes of Harmonics Listed in Figure 25. 
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VI. CONCLUSIONS 


A. FLOW VISUALIZATION 

The original intent of this investigation was to 
produce a flow visualization video of a rotor operating in 
the vortex ring state. The primary motivation was to 
improve upon the video images captured by the classic Drees 
and Hendal smoke video from 1951. [Drees, 1951] Filming the 
flow field on video proved difficult with the limited 
filming and lighting assets. Although, strobe lights 
operating at twice the rotational speed of the two-bladed 
rotor provided superb viewing conditions. With further 
effort, there is significant potential for capturing 
valuable images on video. 

Images captured using a digital camera were extremely 
clear. The rotor wake vortices were captured with 
unexpected clarity. A high speed shutter setting of 1/320 
second captured exceptionally clear images of the helical 
vortex wake. Such images were only seen by the eye when 
using a strobe light to illuminate the flow field. 
Fortunately, the high shutter speed camera captured the 
same images using just flood lighting. 

The use of Particle Image Velocimetry should also be 
considered for use to determine the flow characteristics 
during VRS. Seeding the flow field with small particles and 
illuminating them in a single plane would produce 
interesting results. By tracing the particle motion the 
velocity of the flow field can be mapped, providing 
valuable information about VRS flow. 
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B. VIBRATION ANALYSIS 

Data acquired from thrust measurements were analyzed 
using an FFT analyzer. This revealed surprising correlation 
between the water tunnel model and actual flight data. An 
inexpensive water tunnel experiment using a scaled model 
produced valuable test results. Of primary interest is the 
thrust-induced vibration power spectrum. These frequency 
domain data enabled a vibration analysis and identified the 
presence of increased higher harmonics during VRS flow 
conditions. The data recorded generated interest to 
determine if the same increase in harmonics occurred in a 
full scale aircraft. Analysis of Scheiman data revealed 
this to be true. 

Harmonic analysis of the aerodynamic loading indicates 
that higher harmonic content becomes prevalent during 
operation in the vortex ring state and is a contributor to 
aircraft vibrations. For a two-bladed system, 2P harmonic 
is present in virtually all regions of flight. Higher 
harmonics, such as 4P and 6P are typically only present 
during particular maneuvers encountering blade vortex 
interaction or if the rotor system is poorly designed. 

An accelerometer could provide a warning to the 
cockpit as higher harmonic amplitudes are detected. High 
vibrational components to the thrust response may be 
encountered during various flight conditions but if 
experienced during a typical steep approach into landing it 
may be an indicator of the proximity of the VRS boundary. A 
higher harmonic warning device seems worthy of further 
investigation. 
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